Scaling of polymers in aligned rods 
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We study the behavior of self avoiding polymers in a background of vertically aligned rods that 
are either frozen into random positions or free to move horizontally. We find that in both cases the 
polymer chains are highly elongated, with vertical and horizontal size exponents that differ by a 
factor of 3. Though these results are different than previous predictions, our results are confirmed 
by detailed computer simulations. 

PACS numbers: 



The statistics of a flexible polymer molecule in a back- 
ground of randomly placed vertical rods has been con- 
sideredpreviously 0, Q| ■ The problem was first investi- 
gated [21 in a dynamical context in order to better under- 
stand the effects of entanglements. Later, [3| the equi- 
librium statistics were examined, and an elongation of 
the chain in the vertical direction was predicted. As we 
will show using both scaling and simulations, elongation 
does take place, but is actually much larger than was first 
thought. 

We also analyze how this problem is very closely re- 
lated to a polymer in a background of mobile vertical 
rods. We argue that this would be an interesting system 
to investigate experimentally, perhaps being related to 
liquid crystal polymer mixtures. 

This work has overlap with many other problems of 
interest. For example, vortices in superconductors with 
columnar disorder has been extensively studied 0, 01 and 
is related to the current polymer problem, a major dif- 
ference being that the polymer (or vortex) was stretched 
and threaded the entire vertical dimension having a line 
tension, meaning that it is appropriate to use directed 
polymers in this situation [^. This apparently leads to 
completely different scaling behavior p. 

Our findings are also relevant to the study of mix- 
tures of rods and polymers. It is possible to solubilize 
many rod-shaped particles — for example nanotubes 0, 
or more commonly, virus particles, boehmite rods, and 
cellulose nanocrystals — and rod-polymer mixtures have 
been studied experimentally and theoretically in the con- 
text of demixing and liquid crystalline phase transi- 
tions 11, i, ES 111 El El El. Furthermore, it is now 
possible to possible to fabricate arrays of vertical nan- 
otubes on a substrate which allows for the possibility 
of observing the scaling behavior predicted here (the case 
of quenched disorder) . 

We use a three dimensional cubic lattice model to 
study this problem, where the rods randomly occupy ver- 
tical lines, and the chain avoids itself and the rods. This 
is related to the problem of self avoiding chain (SAW) 
that also is excluded from point defects. In that case, 
an elegant argument of Gates and Ball IT^ cleared up 
a decade of controversy by showing that the quenched 
version of this problem, where defects are frozen, gives 
exactly the same statistics for the SAW as the annealed 



version, where the defects are also mobile. They argued 
that a frozen background of uncorrelated obstacles on 
an infinite lattice could be subdivided into very large re- 
gions. The statistics of obstacles in each region are in- 
dependent of each other. Therefore the statistics of the 
polymer chain can be obtained by doing an average over 
all these regions. The statistical weight given to each 
region gives precisely the same result as an annealed av- 
erage. Their argument trivially extends to the case we 
consider here. 

From the above paragraph, we conclude that for an 
infinite sized lattice, the problem of frozen rods gives 
identical polymer statistics to the case where the rods 
are mobile. Below we will analyze how finite size lat- 
tices alter the above conclusion. But first we analyze the 
annealed (mobile rod) problem. 

Denote the probability distribution of the SAW with 
coordinates ri, r2, . . . , tat, as PsAwl'o^i}- When it is 
placed in a random potential ^(r) at temperature T, 
the annealed average probability distribution becomes 



(1) 



The average depends on the statistics of V. First we 
review the case of a completely uncorrelated V 0, 0| , 
in which case the average becomes an on-site attractive 
term between different monomers. However because the 
SAW excludes all configurations where two monomers sit 
on the same site, this cannot alter the probability dis- 
tribution. Therefore an uncorrelated random potential 
makes no difference to the statistics of an SAW. This 
result has been confirmed by simulations _18J . 

Now we turn to the case of columnar disorder. Here the 
average in eqn ^ now leads to an attractive interaction 
of all monomers that have the same x and y coordinates, 
regardless of their separation in the z direction. In this 
case, we expect different scaling behavior than we would 
observe from an SAW. Intuitively, the chain should con- 
tract in the x-y plane. We now analyze the annealed 
statistical mechanics of this problem using a scaling ar- 
gument. 

A dilute gas of rods of density (per unit area) of p 
and at temperature T gives rise to a (two dimensional) 
pressure p = pT. An SAW of N steps placed in a gas 
of these rods will raise the free energy of the rods by 
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excluding rods from the vicinity of the SAW. This is due 
to the pressure exerted by the rods which means that the 
free energy is raised by pA where A is the cross-sectional 
area of the rods in the x—y plane. Therefore the SAW has 
forces acting on it to decrease A. Working in opposition 
to this is the entropy loss of a polymer chain confined to a 
cylinder of radius Rxy which is 19] ^ N/N^y where N^y 
is the number of monomers corresponding to an SAW of 
dimension R^y ^ N^y, where v « .59 is the excluded 
volume exponent in three dimensions. Minimizing the 
pressure term and tube confinement terms gives: 

R.y - 7Vi/(2+i/'') « N°-^'. (2) 

Note that the exponent 0.27 is much less than in the 
pure SAW case — the pure SAW problem is isotropic and 
so ly is as well. This shows that the scaling of the polymer 
chain in this situation differs greatly from a pure SAW. 
To obtain scaling in the vertical direction R^, we note 
that there are "blobs" of size R^y stacked vertically on 
top of each other. Therefore 
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FIG. 1: Mean monomer separation squared as a function of 
arclength defined in eqn|l] both the total (upper curve) and 
the x-y projection (lower curve) for an SAW of 4096 steps and 
vertical rods with p = 0.32 with mobile rods. The best fit to 
the upper graph is line is 2 x .974. 



Rz - R.yN/Nxy ^ iV3/(2+l/-) ^ (i?^^)3 ^ ^.81^ (3) 

This exponent is thrice the x—y exponent, and is larger 
than v for the pure SAW. Therefore for long chains the 
system becomes highly elongated. The overall radius of 
gyration therefore scales the same way as Rz- 

Note that for long chains, the chain density projected 
on to the x-y plane grows algebraically for large N. 
Therefore for large N we expect complete exclusion of 
rods in this region, which is consistent with our initial 
assumption that this was the case. In addition because 
monomers in far away blobs are uncorrelated, and the 
number of blobs grows with iV, single configurations all 
have very close projected x-y densities which will only 
differ from each other by an inverse power of N. This 
self-averaging property implies that we are also justified 
in regarding the polymer as being in an external poten- 
tial, induced by the rods, and then minimizing the free 
energy to find the scale of this potential. Therefore we 
expect our scaling argument to work very well for large 
iV. 

As argued above, this is identical to the scaling ex- 
pected in the case of a frozen background for a system 
of infinite size. However the case of frozen finite size 
systems is harder to analyze through analytical means 
as the above argument does not apply. In general for 
a finite size system, one does expect a difference be- 
tween the quenched and annealed cases. Therefore we 
performed extensive simulations of this system to verify 
the above predictions and better understand the effects 
of finite sizes. 

Chains of length 32 to 4096 were simulated on a cu- 
bic lattice with different densities of vertical rods. First 
we allowed both rods and chains to move according to 
rules that satisfied detailed balance. The polymers were 
moved according to reptation dynamics (20| and the rods 



performed long range moves between randomly selected 
positions also satisfying detailed balance. We checked 
that our program did indeed give a reasonable excluded 
volume exponent in the case of no rods, « .59 and 
performed a variety of other checks. 

We then analyzed the statistics of our longest chain 
and highest rod density, N — 4096 and p — .32. We plot 
the statistics for the mean monomer separation squared 
as a function of arclength and the same thing but just 
projected in the x-y plane, that is 

(r^(.)).^J,^,.^(|r,-r,n (4) 

These two quantities are shown in fig. the upper and 
lower curve being (r^(s)) and {r'j_{s)') respectively. The 
best fit slope to (r^(s)) is 2 x .974. The scaling theory 
predicts a fractal dimension of 1, which corresponds to a 
slope of 2 on this plot. (r^(s)) has a large plateau indi- 
cating an infinite fractal dimension, consistent with our 
scaling analysis of chains in a tube-like conformations. 

Next we analyzed different chain lengths and rod den- 
sities. The results are shown in figure [21 The x's with 
the upper fitted line going through them represent a fit 
of the exponent for chains of length 128 to 4096 when the 
occupation fraction of rods is p — .32 . The fit gives an 
exponent of .82 which compares well with the predicted 
value .81 given in eqn. (Q. The *'s and D's are the same 
data for filling fractions of .16 and .08 respectively. 

At all filling fractions there is a crossover from small 
chain lengths, which show the usual self avoiding behav- 
ior, to the elongated regime studied here. From the fig- 
ure it is clear that for small chain lengths the radius of 
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FIG. 2: Radius of gyration squared versus chain length on a 
log-log plot for annealed SAW's and vertical rods. The best fit 
line is 2 X .82 for the points at the highest rod filling fractions 
(.32). The -l-'s show the radius of gyration in the x-y plane 
for the same conditions. The other points show the radius of 
gyration at lower rod filling fractions. 



FIG. 3: Radius of gyration squared versus chain length on 
a log-log plot for quenched SAW's and vertical rods. The x 
and square symbols denote a rod filling fraction of .16 and 
.08 respectively. The best fit line is 2 x 1.01 for the points 
at filling fractions of .16. The +'s and *'s show the radius of 
gyration in the x-y plane for the same conditions. 



gyration becomes independent of rod filling fraction as 
expected. It is also apparent from the figure that this 
crossover chain length decreases as rod filling increases, 
which is in agreement with our theoretical interpretation 
of the data, which predicts this cross-over chain length 
Nc scales as Nc p"^"- 

The symbols with the straight line fit going 

through them is the radius of gyration in the x-y plane. 
It is not clear if this has reached an asymptotic power 
law regime even at the larger chain lengths because the 
slope of this line is quite small and therefore crossover 
effects from short chains should be much larger than for 
Rz- The slope of the fit shown is .23 which is nonetheless 
not too far from the predicted value of .27. 

Now we turn to the results found when the rods are 
frozen to random positions. In this case we did not go 
higher than a rod filling fraction .16 because we wanted 
to stay far away from the percolation transition and the 
break down of ergodicity associated with it. The lattice 
size was 128 x 128 in the x-y direction. The chains were 
averaged over a large number of steps. The most time 
consuming case were the longest chains, N = 4096. Be- 
fore averaging, a chain was equilibrated for 1.6 x 10® steps 
and averaged over 1.28 x 10® steps. We then averaged the 
results over 128 realizations of the quenched rods. The 
radius of gyration squared is plotted as a function of the 
number of monomers. The x's and the D's denote a rod 
filling fraction of .16 and .08 respectively. The slope of 
the line going through the .16 filling fraction data gives 
Rg ~ N^-'^^, which is higher than the annealed case. On 
the other hand, the radius of gyration in the x-y plane 
levels off as shown by the two lower curves. The -l-'s and 
the *'s are for a rod filling fraction of 0.16 and 0.08 re- 
spectively. Here, in contrast to the annealed case — figure 
12 — a leveling off in the curves is apparent. 



The higher slope of Rg{N) for quenched as opposed to 
annealed disorder can be understood as follows. As dis- 
cussed above, for large enough lattices one expects the 
quenched average to equal the annealed one. The poly- 
mer migrates around the lattice and finds rare locations 
where it will spend most of its time. In these locations 
there will be a low rod density, and the radius of the de- 
pleted region scales as in eqn. |3 As the chain length 
grows, this implies that these regions must become suc- 
cessively rarer because the size of these holes must grow 
in order to accommodate the chain. On the other hand, if 
the lattice size is finite, then the strength of the most at- 
tractive region is now bounded and cannot grow with N; 
hence the leveling off of the horizontal radius of gyration 
seen in figure |21 A polymer chain in a fixed tube has a 
overall radius of gyration Rg on N 19 ] in good agreement 
with our findings. 

To confirm this explanation, we ran the simulation for 
chains of length 4096 for different horizontal system sizes, 
16 X 16 up to 128 X 128. The horizontal radius of gyration 
for different lattice sizes is plotted in figure0]as a function 
of inverse system size. The point at zero represents the 
annealed cases which should be the same as the infinite 
lattice. As the system width decreases, the horizontal 
radius of gyration increases indicating that the polymer 
is less tightly bound to rod-depleted regions. 

For many chain systems the rods will induce an attrac- 
tive interaction between different chains. The annealed 
case is simplest to understand. When two chains bind, 
they will stack on top of each other but should not inter- 
penetrate as this strongly raises the free energy. From 
the above scaling argument, eqn [51 the free energy of two 
chains binding should scale as the horizontal area occu- 
pied by a chain Rly - j\i2/i2+i/,.) ^ ^o.54^ Therefore 
as polymer concentration is increased from zero, one ex- 
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FIG. 4; Horizontal radius of gyration squared versus horizon- 
tal system size. The horizontal axis is 100 times the inverse 
width. 

pects that there will be a transition between a gas phase 
of isolated chains and a phase of chains stacked on top 
of each other. Eventually it will reach another regime 
where one needs to consider the finite extent of rods in 
the vertical direction. As the concentration is increased, 
at some point different groups of stacked chains will want 
to coalesce. At this point the system will demix into rod 
rich and polymer rich regions. The exact nature of the 
phase diagram is probably quite complex and requires 
a detailed understanding of the behavior of semi-dilute 
polymer solutions in pores |3 1^ j which itself is quite 
dehcate. 

This is also hard to simulate quantitatively (even in 
the annealed case) because the center of mass diffusion 
of polymers is very slow on a monte-carlo time scale. This 
means that for a large system volume, it takes a long time 
for two chains to coalesce. However, it was observed that 
in the pillar regime, two chains that finally managed to 



diffuse into the same tube stayed bound to each other. 

These results found here differ strongly from those 
found for a model for an SAW in the presence of a so- 
lution of vertically aligned short rods where their 
SAW's were of order 100 units but the rods were only one 
or two lattice spacings. In that case, the authors found 
that the polymers went into disk-like configurations in- 
stead of the rod-like ones found here. Here it is worth- 
while examining the case of rods of length L and chains 
whose vertical extent is much larger. In this case eqn. ^ 
implies an attractive potential between monomers with 
the same x-y coordinates whose vertical separation is less 
than L. At separations larger than L we expect that 
chain segments will no longer be stacked in the same tube 
because of the absence of attraction at this length scale. 
This, according to eqn|31 corresponds to a crossover chain 
length of Nc i(2+i/'^)/3 Beyond this separation, differ- 
ent monomers become unbound and we expect that in 
the horizontal direction segments execute a random walk. 
Therefore for large N » N^, R^y L^'^{N/NcY''^ - 
^i/3^_(i+i/(2.))/3^i/2 _ ^-i/(6,.)^i/2^ Therefore we 

expect a slow drift of the SAW in the vertical direction 
that decreases with increasing L. 

In conclusion, we have shown that a polymer in the 
presence of either mobile or frozen vertical rods forms 
tube-like vertical conformations where the tube diameter 
and tube length scale as in eqns [3 and |21 It should be 
possible to observe this unusual behavior experimentally. 
One possible system for the case of frozen rods would 
be in forests of vertical nanotubes Ts'l suitably coated 
to make them repulsive to polymer chains 0. These re- 
sults are also relevant to the more complicated case of 
polymer-rod mixtures, which are usually analyzed assum- 
ing no change in polymer statistics due to the presence 
of rods. 
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